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We study disorder induced transition between quantum Hall states and insulator state in a lattice 
model. We find the positions of the extended states do not change as disorder strength is varied. As 
a consequence, there are direct transitions from all Landau level quantum Hall states to insulator 
states, in contrast to the global phase diagram from early studies based on continuous models. We 
also provide the microscopic understanding of the transition in terms of the topological properties 
of the system. 
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R is well known that all electrons in a two-dimensional 
system are localized in the absence of magnetic field ac- 
cording to the scaling theory of localization [1]. When the 
two-dimensional electron system is subject to a strong 
perpendicular magnetic field, the energy spectrum be- 
comes a series of impurity broadened Landau levels. The 
external magnetic field breaks the time-reversal symme- 
try and as a consequence, extended state appears in the 
center of each Landau band while states away from the 
Landau band centers are still localized. This gives rise 
to the integer quantum Hall effect [2]. Recently there 
have been much interest to understand the evolution of 
the extended states as the magnetic field goes to zero or 
equivalently as the disorder increases such that eventu- 
ally all the extended states should disappear [3-7]. 

Earlier work by Laughlin [8] and Khmelnitskii [9] con- 
cluded that the extended state energy for the nth Landau 
level (LL) behaves in the following form: 

El = {n+ll2)nuj^-±^^, (1) 

where ujc is the cyclotron frequency and r is impurity 
scattering time. According to this formula, E"^ deviates 
from the linear magnetic field dependence when uJcT ~ 1 
and starts to float up as B decreases or as disorder in- 
creases. Based on the floating up picture from contin- 
uous models by Laughlin and Khmelnitskii, Kivelson, 
Lee and Zhang [10] proposed the global phase diagram 
for transitions between the quantum Hall states (QHS) 
and insulator (or localized) states. There are two im- 
portant consequences from the global phase diagram: (i) 
as disorder increases for a fixed field, no direct transi- 
tion between higher LL [n ^ 0) QHS and insulator is 
allowed; (ii) as B decreases for fixed disorder, reentry be- 
havior should be observed in quantized Hall effect, e.g. 

In an earlier work [11] by us and Niu, we found that 
the floating-up picture is not valid in the lattice model. 
We concluded that: (1) the extended state energy Ec 



for each Landau level is always linear in magnetic field; 
(2) for a given Landau level and disorder configuration 
there exists a critical magnetic field Be below which the 
extended state disappears; (3) the lower LLs are more 
robust to the metal-insulator transition with smaller Be- 
We attributed the above results to strong LL coupling 
effect . 

In this paper, we expand our study on disorder induced 
transition between quantum Hall states and insulator 
state in a lattice model. By carrying out finite-size scal- 
ing calculations of localization length, we demonstrate 
that the positions of the extended states do not change 
as disorder strength is varied. As a consequence, we show 
unambiguous evidence from our finite-size scaling analy- 
sis that there are direct transitions from all Landau level 
quantum Hall states to insulator states, in contrast to 
the global phase diagram [10]. We also provide the mi- 
croscopic understanding of the transition in terms of the 
topological properties of the system. 

In the following, we briefly outline our model [12,13] 
and technique. We model our two-dimensional system 
in a very long strip geometry with a finite width (M) 
square lattice with nearest neighbor hopping. Periodic 
boundary condition in the width direction is used to get 
rid of the edge extended states. The disorder potential 
is modeled by the on-site white-noise potential Vim (*' 
denotes the column index, m denotes the chain index) 
ranging from —W/2 to l/F/2. The effect magnetic field 
appears in the complex phase of the hopping term. The 
strength of the magnetic field is characterized by the flux 
per plaquette ((/)) in unit of magnetic flux quanta [4>o = 
hc/e). The Hamiltonian of this system can be written 
as: 

M 

7i = '^'^Vim\im >< im\ (2) 

i m = l 

+ 1*"^ X -^'"1 + ^^irn Jn Ii" >< ' 
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where < im;jn > indicates nearest neighbors on the lat- 
tice. The amplitude of the hopping term is chosen as the 
unit of the energy. For a specific energy E , a transfer 
matrix Ti can be easily set up mapping the wavefunction 
amplitudes at column i — 1 and i to those at column i+l, 



i'i+i 



Ti 



ipi- 



Hi-E -I 
I 



(3) 



where Hi is the Hamiltonian for the ith column, / is a 
M X M unit matrix. Using a standard iteration algo- 
rithm [14], we can calculate the Lyapunov exponents for 
the transfer matrix Ti . The localization length Xm (E) for 
energy E at finite width M is then given by the inverse of 
the smallest Lyapunov exponent. In our numerical cal- 
culation, we choose the sample length to be over 10'* so 
that the self-averaging effect automatically takes care of 
the ensemble statistical fiuctuations. We use the stan- 
dard one-parameter finite-size scaling analysis [15] to ob- 
tain the thermodynamic localization length ^. The scal- 
ing calculations are carried out with varying disorder W 
at fixed energies. According to the one-parameter scal- 
ing theory, the renormalized finite-size localization length 
Xm I M can be expressed in terms of a universal function 
of M/^, i.e.. 



M 



f 



M 



(4) 



where f{x) cx l/x in the thermodynamic limit (M co) 
for localized states while approaching a constant (~ 1) 
when ^ diverges. 
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FIG. 1. Finite-size localization length (Am) as function of 
energy E at magnetic field = 1/8 for disorder W = I 
(circle), W = 3 (triangle), W = 7 (diamond) and W = 12 
(square). The filled symbols are the data for system size 
M = 32 and the open symbols are for M = 64. 



In Fig.l, we present the results of disorder W depen- 
dence of finite-size localization length Xm ■ Fig.l shows 
the Xm at magnetic field cf) = l/H with system sizes 
M = 32 (filled symbols) and M = 64 (open symbols) 
for disorder strengths W = 1 (circle), IF = 3 (triangle), 
W = 7 (diamond) and FF = 12 (square). Due to the sym- 
metry of the lattice model, only the lower energy branch 
results are shown here. In the case of weak disorder with 
FF = 1, we can clearly see that the energy band does 
break up into small subbands (Landau bands). For a 
finite-size system, the states close to the Landau centers 
are extended since their localization lengths Xm depend 
on the system sizes as shown in Fig.l. We have calcu- 
lated Xm for system sizes M = 16, 24, 32, 48, 64 and 
84. We found that Xm at the centers of Landau bands 
scale with the system size M. Therefore, the maxima in 
the finite-size localization length are the locations of the 
extended states even at thermodynamic limit. For those 
states between Landau bands, their localization lengths 
are independent of system size, and are localized states. 
We have done systematic studies with many diffrent FF 
values and found that as FF increases, the peak posi- 
tions in Xm do not shift. Thus, we did not observe the 
fioating-up of the extended state energy for increasing 
the disorder. 

The transport property of a state is determined by the 
topological property-Chern number of the state [16,17]. 
For an eigenstate |m), the boundary phase averaged Hall 
conductance takes the form [16,17] 



' xy/ 



1 

4^ 



92TTy{ei,e2) = C{m)eyh, (5) 



where C{m) is an integer called the Chern number of the 
state |m). States with nonzero Chern numbers are ex- 
tended states which carry current and states with zero 
Chern numbers are the localized states. The transition 
for a state from being extended to localized is caused 
by the cancellation of Chern numbers and is equivalently 
a transition of Chern number from a non-zero integer 
to zero. For the case presented in Fig.l, the extended 
state close to the band center (E = 0) carries Chern 
number C = —3 while all the other extended states 
(E ~ —1.0, —2.0, —3.3 in Fig.l) carry Chern numbers 
C = -|-1 [13,11,18]. As disorder strength FF increases, the 
C = —3 extended state moves toward the band bottom to 
cancel the C = -|-1 extended states one by one. Another 
possibility is for the C = +1 states to move toward band 
center to annihilate the C = —3 state which will cor- 
respond to the fioating-up picture. However, the latter 
possibility is not consistent with our numerical findings. 
It is believed that extended state carrying non-zero Chern 
number can only move from Ei to E2 if there are no lo- 
calized states between Ei and E2, or loosely speaking 
no mobility gap between Ei and £"2. Fig.l clearly con- 
firmed this belief. As disorder increases from FF = 1 to 
FF = 3, all the states close to the band center E = are 
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extended, including those which are localized at W = 1. 
This can be seen by the separation in Xm for M = 32 and 
M = 64. However, those states close to the band bottom 
E = —4 are still showing Landau bands. At = 7, the 
Landau band structure is totally gone and all states for a 
given system size M have similar values for Xm ■ Notice 
that in this case the Xm is of the same order as M, thus, 
all the states are some what extended in nature. Further 
increasing disorder to = 12 will cause all the states 
to be localized and their Xm are much smaller than sys- 
tem size M . Above discussions about the cancellation of 
the Chern numbers is valid for a given finite system. In 
order to understand the transition in an infinite system, 
one has to study the thermodynamic localization length 
^ which is the focus of Fig. 2 in the paper. 



(a) (f)=f/8 




0.0 4.0 8.0 12.0 16.0 
W 




0.0 4.0 8.0 12.0 16.0 
W 

FIG. 2. Thermodynamic localization length { for (a) 
E = —2.0 (diamonds) and —3.3 (down-triangles); and 
{h)E = —0.7 (circles), —1.6 (squares), and —2.8 
(up-triangles). The scaling functions for E = —3.3 (a) and 
E = —2.8 (b) are shown in the insets. 

In Fig. 2 we present our results on thermodynamic lo- 
calization length ^ obtained from finite-size scaling cal- 
culations. To avoid the severe finite-size effect, we have 



used large system sizes in the scaling calculations with 
M = 16,24,32,48,64 and 84. We carry out the scal- 
ing analysis with varying disorder W at fixed energy E. 
In the insets, we show the scaling functions ^'"j^'^ = 
/ (jpTj) for both E = -3.3(a) and E = -2.8(b). The 
scaling plots are quite similar for the other energies we 
have calculated. One can see that all the data fall onto 
one smooth scaling curve from which we obtain the ther- 
modynamic localization length ^. In Fig. 2, we show ^ for 
five different energies E = —0.7 (circles), —1.6 (squares), 
—2.0 (diamonds), —2.8 (up-triangles) and —3.3 (down- 
triangles). 

Presented in Fig. 2(a), E = —3.3 and —2.0 are the 
energies at the centers of the lowest and the first Landau 
bands when the disorder is weak such that the Landau 
bands are well defined as shown by the circles in Fig.l. 
These states are extended with infinite ^ for W < Wc 
and become localized when W > Wc- The important 
conclusion here is that the transitions to the insulator 
states are identical for both Landau bands. Therefore, we 
clearly see the direct transition from high LL quantum 
Hall state to insulator which is not allowed according to 
the global phase diagram [10]. We found similar critical 
disorder Wc (— 4.5) for both Landau bands. However, 
we should mention that it is difficult to estimate Wc for 
the following reasons. At first, the energy value for the 
Landau band center can not be determined very accu- 
rately. Secondly, ^ only approaches infinity from one side 
whereas for the case we are going to discuss below where 
^ diverges from both sides. 

Presented in Fig. 2(b), E = -0.7, -1.6 and -2.8 are 
localized states at weak disorder limit as shown by the 
circles in Fig.l. These energy values are chosen because 
each of them is between two adjacent Landau bands and 
they are the most localized states in the valley. As W 
increases, these states become more "extended" with in- 
creasing ^. At critical disorder strengths Wc, these states 
become truly extended with divergent ^. Beyond Wc, 
they become localized again. The critical strengths Wc 
are different for different energies and it is smaller for the 
energy closer to the band center(£' = 0). Our best esti- 
mations are Wc — 3.5, 4.0 and 4.5 for E = —0.7, —1.6 
and 2.8. As we have mentioned above, the insulator tran- 
sition is caused by the cancellation of the Chern number, 
and Chern number can only move from one energy to 
another if there is no mobility gap in between. The sig- 
nificance of the critical Wc{E) is to allow the negative 
Chern number around band center £" = to move be- 
yond E, thus, to annihilate the -|-1 Chern number state 
on the other side. Once the -|-1 Chern number is anni- 
hilated to zero, that quantum Hall state is transformed 
to an insulator. Therefore, because the difference in Wc 
for different energies, we conclude that less disorder is 
required to induce higher LL to insulator transition. 

Before conclusion, we would like to comment on the 
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validity of the lattice model in describing real experi- 
ments. All the solid state systems should be described 
by lattice models. However, if the relevant length scale 
is much larger than the lattice spacing a, then a contin- 
uous model may be applicable. In a quantum Hall sam- 
ple, the magnetic length I is the relevant length which, in 
general, is much larger than the lattice spacing. On the 
other hand, there is an intrinsical difference between con- 
tinuous and lattice models in this problem in the sense 
that there is no negative Chern number state at finite 
energy in a continuous model. Therefore, most likely, Ec 
will fioat up in a continuous model as described in Eq.(l). 
The natural consequence is the global phase diagram [10] 
which predicts that there exists no direct phase transi- 
tion between high LL quantum Hall state and insulator 
state. Clearly, this type of transition is allowed in the 
lattice model discussed here. In the previous work by 
us and Niu [11] with the same lattice model, we found 
that essential physics of the lattice model prevails even 
at very small field with large l/a. Recent experiment 
[6] might have already observed such transitions between 
high LL quantum Hall states and insulator states. We 
believe that the lattice model is necessary to understand 
this type of experiments for the reason that mentioned 
above. 

After the work was completed, we became aware of 
a recent work by Yang and Bhatt [19]. They obtained 
similar numerical results as ours, however, their interpre- 
tation of the results are quite different. They conclude 
that although the positions of the extended states do not 
change much as were found in our work, the band posi- 
tion moves down as disorder W increases. Therefore, rel- 
ative to the band bottom, the positions for the extended 
states fioat up. There are two problems with this argu- 
ment: (i) There have been calculations [11] with changing 
magnetic field at fixed disorder. In that case, band bot- 
tom does not change and, even there, the fioating-up of 
Ec is not seen, (ii) According to the original fioating-up 
picture put forward by Laughlin and Khmelnitskii, the 
fioating up only happens when disorder W is of the order 
of Landau level spacing. However, using the reasonings 
by Yang and Bhatt, one would say that there is fioating- 
up even for infinitesimal increasing in disorder W which 
is quite different from the conventional understanding of 
the fioating-up picture. 

In summary, by calculating numerically the localiza- 
tion length in a lattice model, we studied the disorder 
induced transition between quantum Hall state and in- 
sulator state. We find the positions of the extended states 
do not change with varying disorder strength. As a conse- 
quence, there is direct transition from high Landau level 
quantum Hall state to insulator state, in contrast to the 
global phase diagram based on continuous models. This 
transition can be understood from the cancellation of the 
topological Chern numbers. 
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